arXiv:1508.00564vl [hep-ph] 3 Aug 2015 


IPPP/15/48, DCPT/15/96, CP3-15-23, MCnet-15-20 


Dark matter production through loop-induced processes at the 
LHC: the s-channel mediator case 

Olivier Mattelaer 1 , Eleni Vryonidou 2 

1 Institute for Particle Physics Phenomenology (IPPP), Durham University, Durham, DH1 3LF, 

United Kingdom 

2 Centre for Cosmology, Particle Physics and Phenomenology (CP3), Universite catholique de Louvain, 

B-1348 Louvain-la-Neuve, Belgium 


Abstract. We show how studies relevant for mono-X searches at the LHC in simplified models featuring a 
dark matter candidate and an s-channel mediator can be performed within the MadGraph5_aMC@NLO 
framework. We focus on gluon-initiated loop-induced processes, mostly relevant to the case where the 
mediator couples preferentially to third generation quarks and in particular to the top quark. Our imple¬ 
mentation allows us to study signatures at hadron colliders involving missing transverse energy plus jets 
or plus neutral bosons ( 7 , Z, H), possibly including the effects of extra radiation by multi-parton merging 
and matching to the parton shower. 


1 Introduction 

Searching for Dark Matter (DM) candidates is a leading 
priority of Run II at the LHC. As the DM particles are ex¬ 
pected to be stable and neutral under the Standard Model 
(SM) interactions, they do not leave any trace in the de¬ 
tectors. Therefore, the DM searches focus on the associ¬ 
ated production with visible SM particles. These processes 
lead to missing transverse energy and are often classified 
as mono-X searches, where X can be jet(s), photon, Z, W 
or even a Higgs boson. 

DM model building is a very rich and exciting theo¬ 
retical activity, with many mechanisms and models be¬ 
ing explored in the community. Such a wealth of options, 
however, makes top-down searches somewhat limited in 
scope and calls for more general bottom-up approaches. 
These are generally classified into two groups. The first 
involves effective operators, describing the direct inter¬ 
action between the DM particles and the SM ones, via 
higher dimensional operators. These so-called Effective 
Field Theory (EFT) models arise in the presence of suffi¬ 
ciently heavy additional states, not directly accessible in 
our experiments, that can be integrated out leading to a 
contact interaction between the dark matter and the SM 
particles. In this case, experimental searches can set limits 
on the scale A (assuming unit coupling) which determines 
the strength of the interaction. 

The second class of models, called simplified models of 
DM, is based on a single particle exchange between the 
DM and SM particles. The particle exchanged is called 
the mediator, a neutral boson in the case of s-channel 
mediator, or a more generic state in the case of f-channel 
ones. Simplified models generally include more parameters 
than EFT models, as the couplings of the mediator to the 


DM and SM particles, the width and mass of the media¬ 
tor play a role in the phenomenology, and hence they can 
parametrise more reliably the kinematics of the produc¬ 
tion process when the mediator is not sufficiently heavy 
to be integrated out. 

Both classes have been extensively studied in the liter¬ 
ature. Simplified models for DM are presented for example 
in PQ, while a review of simplified DM models and their 
relevance for missing energy searches at the LHC is given 
in 0. The EFT limit of these interactions for the scalar 
mediator case was studied recently in [2, where limits on 
the EFT operators are set, based on the missing transverse 
energy searches at the LHC. Various other works focus on 
DM production in the EFT, providing results at NLO in 
QCD for various mono-X processes BH3- We refer the in¬ 
terested reader to the recent reviews in Refs. 00 for a 
more complete list of models and processes to be studied 
at the LHC. 

Within both classes of models, an especially interest¬ 
ing scenario arises when either the dark matter or the 
mediator only couples to the third generation quarks and 
in particular to the top quark. In this case, the only DM 
production mechanism at tree-level becomes production 
in association with heavy quarks, while all other processes 
have to proceed through heavy quark loops in gluon fu¬ 
sion. Due to the high gluon luminosity at the LHC, these 
processes can be important and have already received at¬ 
tention in the literature. Within the EFT models, this 
scenario has been studied in 0HU] for the mono-jet sig¬ 
nal for scalar or pseudoscalar contact interactions. In the 
same EFT setup, m considers the production of DM in 
association with two jets, whose angular correlations can 
provide a handle to distinguish between scalar and pseu¬ 
doscalar interactions. In the context of a simplified model, 
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the mono-jet process from top quark loops for a scalar and 
pseudoscalar mediator has also been discussed [T2[[T3] . 

In this work, we present results of the implementation 
of a simplified DM model in the MadGraph5_aMC@NLO 
framework (MG5 aMC henceforth) |T|], including a Dirac 
fermion as dark matter and either a scalar, pseudoscalar, 
vector or axial-vector mediator. In particular, we focus on 
the case where the mediator couples only to the top quark 
(and to bottom in the case of an axial vector to avoid 
anomalies), and study the loop-induced contributions to 
the most commonly searched for mono-X processes for 
various masses of the mediator and the dark matter par¬ 
ticles. 

This paper is organised as follows. In section [2] we 
briefly describe the simplified model employed, and in sec- 
tion[3]the calculation setup. In section[l]we present results 
for the jets plus missing transverse energy signature com¬ 
paring the full top-quark mass dependence with the top- 
EFT approach while in section [5] we discuss mono-Z, Higgs 
and photon processes. We draw our conclusions in section 

EE 

2 Simplified Model 

We assume a simplified model where the mediator only 
couples to top quarks and to the DM particles. The DM 
particle (x) is taken to be a Dirac fermion in this study 
but the implementation of the model (see next Section) 
is flexible to allow a real or complex scalar. The mediator 
can be chosen to be either a scalar (Y 0 ) or a vector (Yj). 
In the scalar case, the interaction Lagrangian is given by: 

^dm ~ x(9dm + *5dm7 5 )x^o, 

+ *5t P 7 5 )i Y 0 . (1) 

For convenience we normalise the scalar and pseudoscalar 
interaction with the top quark (t) by the top Yukawa cou¬ 
pling (y 4 = 2^). In the vector case, the interaction La¬ 
grangian is given by: 

£dm = X7 ij.(9dm + 9dm1 5 )x^i > 

C Vm = tibaY + 9?l 5 )tY? + bn{~9tl 5 ) bY i ■ ( 2 ) 

For the axial-vector case, we also introduce a coupling to 
bottom quarks (61 and fixed it to be opposite to the one of 
the top quark (gf) to ensure the cancellation of the gauge 
anomaly. 

The couplings, mass and width of the mediator and 
the mass of the DM can be extracted from the predictions 
of more concrete DM models. In particular, the width of 
the mediator can be computed by considering the particle 
content of Csm + £ dm and possibly of a more involved 
dark sector. 

3 Method and tools 

3.1 Technical setup 

The computation is performed within the MG5 aMC frame¬ 
work. For this study we employ the latest version (2.3.0), 


which allows automatic event generation for loop induced 
processes [15]. Within MG5aMC, MadLoop m com¬ 
putes the one-loop amplitudes, using the OPP integrand 
reduction method m (as implemented in CutTools USD- 
We use the UFO model of m , generated by FeynRules/ 
NLOCT I201I2T] . The model can be downloaded from the 
FeynRules catalogue of models [22]. We stress here that 
the implementation of the model has been validated by 
comparing to the SM. This was achieved by setting various 
parameters to the corresponding SM values and studying 
processes with a Higgs and a Z , mimicked by a Yo and Yi 
respectively, and perfect agreement has been found. 

As the processes we consider in this work start at one 
loop (at leading order), NLO corrections are not avail¬ 
able. Such NLO computation would require some chal¬ 
lenging two-loop multi-scale integrals, most of which are 
not yet available. In order to provide a better description 
of the kinematics, one can employ the method of Matrix- 
Element- Parton Shower (ME+PS) matching/merging [33] . 
ME+PS schemes allow the consistent combination of ma¬ 
trix elements with different jet multiplicities via their match¬ 
ing to a parton shower. Merging of samples of different 
multiplicities in loop-induced processes has been done within 
MG5 aMC in [23] for H+jets and more recently in [21] 
for ZH associated production, by employing a reweight¬ 
ing procedure. ME+PS results for Higgs plus jets, ob¬ 
tained automatically with the interface of MG5 aMC to 
Pythia6 [26] have been presented in [15]. 

The implementation of ME+PS method in MG5 aMC 
comes in two variants: the traditional fcy-MLM and the 
shower-fey schemes. The two give comparable results as 
discussed in 122 ■ In this study we will employ the tra¬ 
ditional fe T -MLM scheme, and in particular the most re¬ 
cent implementation of the scheme in conjunction with 
Pythia 8 j28[]29], for events generated in MG5 aMC. 

The value of the merging scale Q cut is selected on a 
process-by-process basis to ensure that there is a smooth 
transition between the ME and PS regimes. In practice, 
this is assessed by examining the differential jet rate dis¬ 
tributions, which show whether the transition is indeed 
smooth. The differential jet rates are obtained by cus¬ 
tomising the Pythia8 routines. The other distributions 
are obtained by passing the merged samples through Mad- 
Analysis5 HSU, which is interfaced to Fast Jet [22 
for jet reconstruction. For the jet clustering we use a min¬ 
imum jet transverse momentum of 25 GeV, and employ 
the anti-fey algorithm [33] with a radius of R = 0.4. 

While it is straightforward to study ME+PS merging 
for all mono-X processes considered here, for the sake of 
brevity and simplicity we will only present merged sam¬ 
ple results for the jets and missing transverse energy sig¬ 
nal, while for the rest of the processes we will only show 
parton-level results. 

For the results presented here, we use the MSTW2008LO 
]331 parton distribution functions (PDFs) and the cen¬ 
tral renormalisation and factorisation scales are set to 
half the sum of the transverse mas ses of the final state 
particles: p° = fi°R = = \ \j m 'i + Pt %■ I 11 our re ~ 

suits, scale variations are obtained by varying the scales 
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independently in the range of p°/2 < hr,f < 2/r°, com¬ 
puted automatically with SysCalc [35] . All results will 
be presented for the LHC with 13 TeV centre-of-mass en¬ 
ergy. The computation is performed within the 5-flavour 
scheme, unless otherwise stated. 

3.2 Benchmark points 

In order to demonstrate various physics effects within this 
model, we concentrate on a series of benchmarks. In all 
cases, we assume the same nature of coupling between 
mediator-DM and mediator-tops and concentrate on four 
scenarios: the scalar one (gf, M = gf — 1) , the pseudo¬ 
scalar one (g£ M = g[ = 1), the vector one (gY M = gY = 
1) and the axial-vector one (g^ M = gY = !)■ The non- 
specified couplings in these four scenarios are set to zero. 
We also select three mass benchmarks, shown in Table Q] 
We find this choice illustrative as these represent three 
different physics cases. The first benchmark dubbed ‘ res¬ 
onant' corresponds to the resonant production of the me¬ 
diator and its decay into a pair of dark-matter particles. 
The second benchmark ‘heavy med^ator , involves a heavy 
mediator and a light dark matter particle. The width asso¬ 
ciated to the mediator is, in that case, quite large and in¬ 
teresting off-shell effects can be observed. The final bench¬ 
mark ‘ heavy DM ’ corresponds to the case where the me¬ 
diator can not be produced on-shell, as its mass lies below 
the 2 m x threshold. 

For completeness, we show the mediator widths ob¬ 
tained for the three benchmark points and for the four 
possible couplings scalar (S), pseudoscalar (P), vector (V) 
and axial-vector (A) in Table [5] In the first scenario the 
mediator decays only to DM, in the second to DM and 
top quarks and in the third to top quarks only. For the 
axial vector mediator decay into bottom quarks is allowed 
and included in the width computation for all benchmarks. 
The computation of the width has been performed at tree- 
level with MadWidth [5B1 . 

4 Jets and missing transverse energy signal 

4.1 Total cross-section results for the mono-jet process 

The first process we consider is the production in asso¬ 
ciation with QCD jets, which leads to a missing energy 
signature though the production of the mediator and its 
decay to DM particles. A sample of the Feynman diagrams 
contributing to this signal are shown in Fig. [T] 

In this study we will be focussing on the shapes of the 
differential distributions, and therefore we will typically 
show normalised ones. Nevertheless, we find it instruc¬ 
tive to start by presenting the cross-sections we obtain 
for the mono-jet process for the various benchmarks dis¬ 
cussed above, for two cuts on the jet transverse momen¬ 
tum, a minimal cut of 50 GeV and a more realistic one 
of 200 GeV. In this simple process at the parton-level, 
the jet transverse momentum is identical to the missing 
transverse energy, as the DM pair is recoiling against the 


parton/jet. Results are shown in Tabic [3] for the scalar and 
pseudoscalar mediators, including the scale and PDF un¬ 
certainties for the LHC at 13 TeVQ The scale uncertain¬ 
ties can exceed 50% in some cases, while the PDF ones 
rise from ~1% in the resonant case to 2% in the heavy 
DM scenario. Such large scale uncertainties are similar 
to those observed for loop-induced processes within the 
SM , Tf>] . The heavy DM scenario gives the smallest cross- 
sections, as these are suppressed by the DM pair produc¬ 
tion threshold which lies at 1 TeV, leading to sub-fb cross- 
sections. The resonant production is on the other hand 
enhanced, as expected. The gluon fusion amplitudes differ 
for a scalar and a pseudoscalar mediator, by a factor of 
2/3 in the infinite top mass limit. This manifests also here 
with the pseudoscalar mediator giving consistently larger 
cross-sections. We note here that the heavy mediator sce¬ 
nario gives larger cross-sections than the heavy DM one, 
as a significant fraction of the cross-section comes from 
the off-shell region, well below the mass of the mediator 
(1 TeV). These results have been obtained for the cou¬ 
plings set arbitrarily to 1. By modifying these couplings 
the cross-sections can be enhanced/reduced accordingly. 

The corresponding results for the mono-jet cross-sections 
for the vector and axial-vector mediators along with the 
scale/PDF uncertainties for the three benchmarks are shown 
in Table Q] Similar observations regarding the scale and 
PDF uncertainties can be made in this case. The cross- 
sections for the three benchmarks with the vector me¬ 
diator follow the same patterns as those for the scalar 
and pseudoscalar. We notice however that the axial-vector 
cross-sections, especially for the heavy DM scenario, are 
significantly larger than the corresponding vector ones. 
This effect, which will also be evident in the other mono- 
X processes studied in the next section, originates from the 
non-conservation of the axial vector current. The media¬ 
tor propagator contains a term proportional to p^p^/M^ 
that when contracted with the axial-vector current leads 
to terms proportional to mtm x /MY 1 , and leads to an en¬ 
hancement of the cross section for large DM masses. 

4.2 Differential results for the scalar and pseudoscalar 
mediators 

To provide an accurate description of the kinematics be¬ 
yond LO, we consider merged samples of 0, 1 and 2-jet 
multiplicities. For comparison purposes, in the plots shown 
below we also include the main background, Z(—> vv)+ jets, 
obtained for consistency with LO merging of 0, 1 and 2- 
jet samples. Both computations are performed in the setup 
described in the previous section. 

We will compare results for the three scenarios pre¬ 
sented above, showing the distributions for the most rel¬ 
evant kinematical variables. As an example of the checks 
performed to ensure the merging scale chosen was an ap¬ 
propriate one, we show for one process -the resonant scalar 

1 We note here that the corresponding LO results for Z(—> 
uu)+jet are 711 and 11.7 pb for a cut of 50 and 200 GeV 
respectively. 
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Benchmark 

Resonant 

Heavy mediator 

Heavy DM 

Mediator mass 

200 

1000 

400 

Dark matter mass 

50 

1 

500 


Table 1 . Mass benchmarks in GeV. 


Benchmark 

S 

P 

V 

A 

Resonant 

5.17 

6.89 

5.17 

19.3 

Heavy Mediator 

88.0 

94.5 

105.7 

172 

Heavy DM 

3.10 

11.89 

22.2 

36.0 


Table 2. LO widths in GeV for the various mediators. The computation of the width has been performed with MadWidth [36l . 



Fig. 1. Feynman diagrams contributing to jets plus missing transverse energy signal in the simplified model. 


production- the differential jet rate distributions in figure 
[21 This distributions have been obtained with Pythia8 
with a matching scale (Q cut ) of 60 GeV. 

The normalised distributions for the hardest, second 
hardest jet and the missing transverse momentum for the 
three scenarios are shown in Figs. [3] and 0] respectively. 
The invariant mass distributions of the DM pair for the 
various benchmarks are presented in Fig. [5] The resonant 
curve displays a sharp resonant peak, while the heavy- 
mediator scenario has an important tail at low invariant 
masses, with the off-shell region contributing significantly 
to the cross-section. In both curves a threshold effect can 
be observed at 2 mt when the top quarks running in the 
loop become on-shell. We note that the mass of 1 TeV for 
the mediator is not sufficiently high for the EFT approach 
to be valid as the mass of the mediator is probed, as shown 
clearly in Fig. [5j For the third and final scenario heavy 
DM- the production threshold lies above the mediator 
mass and therefore no resonant structure arises. 

We see that the Z+ jets background falls faster than 
any of the DM scenarios, for all the transverse momen¬ 
tum distributions shown. This implies that while inclu¬ 
sively it is overwhelming, deviations from the background 
can be observed more easily by searching in the boosted 
regions. Comparing the three scenarios, we notice that the 
distributions for the resonant scenario, fall more rapidly, 
while the heavy DM and heavy mediator scenarios lead to 
harder distributions in the tails. 

The corresponding results for the pseudoscalar media¬ 
tor are shown in Figs. [6] [7] and [5] These are almost iden¬ 
tical to the results of the scalar propagator, with no vis¬ 
ible difference in the normalised distributions for the jet 
and missing transverse energy. We recall here that a vari¬ 


able that can be used to distinguish between scalar and 
pseudoscalar is the azimuthal separation between the two 
leading jets produced in this process, as discussed for the 
SM Higgs in H + 2 jets in [37] and for DM searches in [11] . 
The slightly sharper threshold visible at 2m* is due to the 
different behaviour of the scalar and pseudoscalar ampli¬ 
tudes. 

For the case of scalar or pseudoscalar mediators, one 
could also consider production cross-sections in the infinite 
top mass limit, i.e. employing the Lagrangian of the form 

m-- 

c = + ^gfG^Yo. (3) 

Such an approach has the clear advantage of being much 
simpler than computing full loop amplitudes and thanks 
to this simplicity, the possibility of including NLO correc¬ 
tions in QCD, see for instance m ■ Nevertheless, the accu¬ 
racy of the infinite top approximation with respect to the 
exact loop computation needs to be assessed on a case-by- 
case basis. As an example, we show, in Fig. [9] a comparison 
for the missing transverse energy distribution in the case 
of the scalar mediator. It is clear from the plot that inte¬ 
grating out the top quark leads to harder distributions in 
the tails and the top-EFT result overshoots the loop-one 
for all three scenarios. We see that for the resonant case, 
the infinite top mass limit provides a reliable prediction 
of the distribution shape up to 200 GeV. These observa¬ 
tions are qualitatively consistent with the corresponding 
studies for Higgs production in the SM, where the infinite 
top mass limit fails at high Higgs transverse momentum. 
Considering the fact that the DM searches focus on the 
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Benchmark 

Scalar 

Pseudoscalar 

Jet pr cut in GeV 

50 

200 

50 

200 

Resonant 

a ■in +49% +0.8% 

^±.11 31 % o.9% 

n o/i/i +50% +1.1% 
u.z^ _ 31% _ 1>1% 

-| rv -i +49% +0.8% 

1U - 1 -31%-0.9% 

n ^Q/i +50% +1.1% 
U.004: _ 31% _ 11% 

Heavy mediator 

q 99 in -2 +55% +1.7% 
O.ZZ • IU 33% -1.7% 

4.92 

qn -3 +55% +1.8% 

' 1U -33% -1.8% 

A OQ in -2 +55% +1.6% 
^.ZO ■ IU —33% -1.6% 

6.47 

in -3 +54% +1.8% 

1U -33% -1.8% 

Heavy DM 

A QQ 1 n— 3 +55% +2.2% 
4.00 • IU —33% -2.1% 

8.54 

in -6 +55% +2.3% 

' 1U -33% -2.2% 

-i iyq -i r\—4 +54% +2.0% 

10 ' 1U -33% -2.0% 

3.35 

in -5 +54% +2.1% 

1U -33% -2.1% 


Table 3. Cross-sections in pb for pp — > XXJ at 13 TeV for the scalar and pseudoscalar mediators for a cut of 50 and 200 GeV 
on the jet transverse momentum and the corresponding scale and PDF uncertainties for the LHC at 13 TeV. 


Benchmark 

Vector 

Axial-vector 

Jet pr cut in GeV 

50 

200 

50 

200 

Resonant 

n /l&T +5i% +1.1% 

U.^Of 3i% -1.2% 

Pi 10/1 +51% +1.4% 

U.1LF4 _ 32% _ 1A% 

i i r +50% +0.7% 

11,0 -31% -0.9% 

n9 +50% +1.1% 

±.uz _ 31% _ 11% 

Heavy mediator 

2.68 

in -4 +56% +2.1% 

‘ 1U -34% -2.1% 

1.55 

in -4 +57% +2.4% 

' 1U -34% -2.4% 

r: ^1 in -3 +52% +0.9% 

0.01 • IU _ 32% _ 10% 

8 97 • 1 0 -4 + 53< ^ + 1 - 4 % 

° 1U -32% -1.4% 

Heavy DM 

1.48 

n -6 +57% +2.9% 

' 1U -34% -2.8% 

1.09 

in —6 +57% +3.0% 

' 1U -34% -2.9% 

19 o in -3+54% +2.0% 
l.ZO • IU 33% —2.0% 

9 ^0 1 0 — ^ +54% +2.1% 
Z.OU • IU _ 33% _ 2 .i% 


Table 4. Cross-section in pb for pp —>■ XXJ f° r the vector and axial-vector mediators for a cut of 50 and 200 GeV on the jet 
transverse momentum and the corresponding scale and PDF uncertainties for the LHC at 13 TeV. 




Fig. 2. Differential jet rate distributions for a scalar mediator 
in the resonant DM production scenario, for a merging scale of 
60 GeV. 


boosted regions to overcome the large SM backgrounds, to 
avoid overestimating the signal and consequently setting 
inaccurate limits on the various DM model parameters, 
one needs to resort to the loop computation. 




Fig. 3. Hardest and second hardest jet transverse momentum 
distribution for pp —> XX+ jets for a scalar mediator. 

4.3 Differential results for the vector and axial-vector 
mediators 

Following the same procedure, one can obtain results for 
the vector and axial vector mediators. When merging sam¬ 
ples of different multiplicities, one needs to take into ac¬ 
count that for a vector mediator, the 0-jet contribution 
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Pt‘“ [GeV] 


Fig. 4. Missing transverse momentum distribution for pp 
XX+ jets for a scalar mediator. 



m xx [ GeV l 


Fig. 5. Invariant mass distribution for the DM pair pp — >• XX+ 
jets for a scalar mediator. 


is zero due to Landau-Yang theorem mm- The quark- 
gluon initiated contributions to the 1-jet process also van¬ 
ish. In fact, any diagram involving a triangle with two 
gluons and a vector mediator as the external legs van¬ 
ishes, due to charge conjugation invariance. Therefore, in 
the vector case, we only merge the 1 and 2-jet contribu¬ 
tions. For numerical stability we have removed the trian¬ 
gle diagrams at the time of generation of the code. The 
differential distributions are shown in Figs. Util fill and [T21 

The three scenarios display the same pattern as for the 
scalar and pseudoscalar, i.e. the resonant case giving the 
most rapidly falling distributions. The difference in the 
shape between the vector and scalar/pseudoscalar distri¬ 
butions originates from the absence of the 0-jet sample. 
Consequently, both the missing pp and hardest jet dis¬ 
tributions go to zero at low pp and present a maximum 
at rather large pp values. The distribution of the invari¬ 
ant mass of the DM pair for the heavy mediator scenario 
shows a peak at both the 2 nit threshold and the mass of 
the mediator. 

Similar results can be obtained for the axial-vector me¬ 
diator. In this case though, there is a subtlety related to 
the fact that if the axial-vector mediator couples only to 
the top quark, this will lead to a gauge anomaly. In the 
SM the gauge anomaly for the Z is exactly cancelled, as 
the axial vector coupling to the Z for up and down quarks 




Fig. 6. Hardest and second hardest jet transverse momentum 
distribution for pp — > XX+ jets for a pseudoscalar mediator. 


differs by a minus sign. Only coupling the mediator to the 
top, would imply that the theory is anomalous, which is 
not expected in any UV complete theory. As we have al¬ 
ready discussed in the introduction, a minimal solution to 
this problem is to allow in this particular case the axial- 
vector mediator to couple also to the bottom quark, with 
a coupling opposite in sign to that of the top. In this case, 
we perform the computation in the 4F scheme to again 
focus on the loop-induced production. 

The differential distributions are shown in Figs. [13] ITT! 
and m The relative shapes for the three scenarios fol¬ 
low the same patterns as the other three couplings for the 
transverse momentum distributions, i.e. the resonant one 
falls more rapidly. A difference compared to the vector me¬ 
diator is seen in the DM invariant mass distribution for 
the heavy mediator scenario. Due to the presence of the 
bottom quarks in the loops, a significant fraction of the 
cross-section lies at low invariant masses. The fact that the 
width of the axial-vector mediator is significantly larger, 
further enhances the off-shell contribution. In fact, both 
the transverse momentum and missing energy distribu¬ 
tions are softer than the corresponding ones for the vector 
mediator. Nevertheless, these distributions remain harder 
than those of the Z background, facilitating boosted tech¬ 
nique based searches. 
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Fig. 7. Missing transverse momentum distribution for pp 
XX+ jets for a pseudoscalar mediator. 



Fig. 8. Invariant mass distribution for the DM pair pp — >• XX+ 
jets for a pseudoscalar mediator. 



Fig. 9. Missing transverse momentum distribution for pp —>• 
XX+ jets for a scalar mediator, using the exact loops (solid 
lines) and the top EFT (dashed lines) and the corresponding 
ratio. 

5 Other mono-X signals 

In this section, we consider the possibility of other mono- 
X signals within our simplified model. Three cases are 
considered: mono -Z, mono-Higgs and mono-photon. In all 
cases the production mode is through top quark loops in 



Fig. 10. Hardest and second hardest jet transverse momen¬ 
tum distribution for pp —> XX+ j e ^ s f° r a vec t° r mediator. 

gluon fusion. The only contribution is via box diagrams 
as the mediators only couple to the top quark. Sample 
diagrams are shown in Fig. [TB] We study the same three 
benchmark points discussed in the previous section. How¬ 
ever, in this section, we limit ourself for simplicity to the 
parton level, and focus on some general physics consider¬ 
ations. 

5.1 mono-Z 

We start by considering the mono-Z associated produc¬ 
tion. The total cross-sections for pp —> XX% f° r unit cou¬ 
plings at 13 TeV are given in Table [5j We find that the 
scalar and pseudoscalar mediators give results of the same 
order of magnitude. Note that, for both scalar and pseu¬ 
doscalar mediator production, it is only the axial-vector 
coupling of the top to the Z boson which contributes, due 
to charge conjugation invariance. 

Similarly to the multi-jet case, the axial-vector medi¬ 
ator gives much larger cross-sections than the vector one. 
As discussed in Section l4~T1 this is due to the propagator 
term p^p 1 ' /My which, when contracted with the axial- 
vector current, leads to terms proportional to m t m x / My 
potentially divergent in the limit of My 1 —> 0, due to the 
non-conservation of the axial current. 

The normalised parton-level missing transverse mo¬ 
mentum distributions, i.e. the transverse momentum of 
the mediator, for the four types of mediator couplings are 
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Benchmark 

scalar 

pseudoscalar 

vector 

axial-vector 

Resonant 

o QQ in -2+36% +1.3% 

• iu _ 25% _ 14% 

q 90 in -2 +36% +1.3% 
o.zo • iu _ 25% _ 14% 

q 9 ft in -3+35%+1.1% 
o.zu • iu _ 24% _ 13% 

OQO in -2+35% +1.1% 

• iu _ 24% _ 12 % 

Heavy mediator 

9 9 /-J I —4 +43% +2.5% 
Z.ZU • 1U _ 28 % _ 2 . 5 % 

9 rjo 1 r»—4 +43% +2.6% 
Z.UO • IU _ 2g% _ 2 . 5 % 

9 1 r 1 r» — 6 +42% +2.2% 
Z.IO-IU -28% -2.3% 

1 1 n ~ 4 + 41% + 2 -o% 

A *°Z ■ 1U -27% - 2 . 0 % 

Heavy DM 

A 7^ 1 n —7+45% +3.5% 
-29% -3.4% 

1 Aft in -6 +44% +3.2% 

1 .^U- 1 U 28% - 3 . 1 % 

1 fK in- 8 + 44% + 3 - 1% 
l.UO • IU 28% 3 0 % 

1 in in- 5 + 44% + 3 - 3% 

1 . 1 U- 1 U 28% - 3 . 2 % 


Table 5. Cross sections (in pb) for gluon induced mono-Z production at the LHC at y/s = 13 TeV for three mass benchmarks. 
A technical cut of 2 GeV has been set on the transverse momentum of all final heavy states. No Z branching ratios are included. 



Fig. 11. Missing transverse momentum distribution for pp —> 
XX+ jets for a vector mediator. 



m xx [ GeV ] 


Fig. 12. Invariant mass distribution for the DM pair pp —>• 
XX+ jets for a vector mediator. 


shown in Figs. [T7] and [TU Again the resonant scenario 
gives a distribution which falls more rapidly, and we find 
that the shape of the distributions are the same for the 
scalar and pseudoscalar mediators. These shapes are simi¬ 
lar to those given by the box contribution to the gg —> HZ 
production is the SM, as shown in m- In the SM though, 
there is a significant, nearly exact at high p-r, cancella¬ 
tion between the box and triangle diagrams, which leads 
to suppressed tails. This cancellation is absent in the DM 
scenarios, and the distributions fall very slowly. We note 
that the total cross-section for the resonant case is smaller 
than the SM one, even though the masses {mu ~ 2 m x ) 
are similar, as it is the triangle contribution that is the 
dominant one in the SM at low energies. 

Distributions corresponding to vector and axial-vector 
mediators differ not only in the normalisation but also in 




Fig. 13. Hardest and second hardest jet transverse momen¬ 
tum distribution for pp —> XX+ jets for an axial-vector media¬ 
tor. 


the shapes for the resonant and heavy DM scenarios. From 
charge conjugation invariance we know that there is no 
mixed term c a c v ( c a /c v are respectively the axial/vector 
couplings of the Z boson), and only the vector and axial- 
vector Z couplings contribute for the vector and axial me¬ 
diators respectively. As discussed in the original computa¬ 
tion of gg —> ZZ in m the contributions from (? v and 
are different and lead to enhancements for massive quarks 
due to the non-conservation of the axial current. 

For this kind of process, it is important to investigate 
the importance of additional QCD radiation, similarly to 
what we observe in the SM for gg —>■ ZH EU- We illustrate 
this effect in Fig.[l!J]by showing in the ‘resonant’ case (200 
GeV scalar mediator) the missing transverse momentum 
for the zero and one-jet multiplicities. We see that these 
contributions are not suppressed compared to the 0-jet 
ones. In the SM this effect, i.e. the importance of the 1-jet 
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Fig. 14. Missing transverse momentum distribution for pp 
XX+ jets for an axial-vector mediator. 



Fig. 15. Invariant mass distribution for the DM pair pp —► 
XX + jets for an axial-vector mediator. 



Fig. 16. Representative Feynman diagrams contributing to 
the mono-X processes in the simplified model. 



Fig. 17. Missing pr distribution for pp — » XX% f° r the scalar 
(S) and pseudoscalar (P) mediators. 



Fig. 18. Missing pr distribution for pp — » XX Z f° r the vector 
(V) and axial-vector (A) mediators. 

contribution, is even more pronounced at high pp as the 0- 
jet amplitude is extremely suppressed by the cancellation 
between box and triangle diagrams. In any case, it might 
be important experimentally to perform inclusive searches 
rather than applying a jet veto, as such a cut will elimi¬ 
nate these enhanced 1-jet contributions. A more accurate 
prediction of the shapes can be provided by LO merging 
and matching to the PS, a procedure that is automatic in 
MG5 aMC. 



Pt‘“ [GeV] 


Fig. 19. Missing pr distribution for pp — » Vo Z(j) for the scalar 
mediator of 200 GeV. The pp —» YoZj distributions are shown 
for two cuts on the transverse momentum of the jet. 


5.2 mono-Higgs 

Studies of mono-Higgs production related to dark matter 
searches, have already been presented in the literature. In 
the scenario studied in msa, the dark matter particles 
only couple to the Higgs and therefore the only production 
mode is in association with a Higgs boson. In addition to 
the mono-Higgs production through the coupling of the 
mediator to the Higgs, the gluon fusion contribution to 
the mono-Higgs production in the infinite top-mass limit 
has been studied in [I]. 

The gluon fusion contribution with the exact top-mass 
dependence can be obtained with our implementation of 
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the simplified model. Mono-Higgs production in this sim¬ 
plified model is not possible with a vector mediator, sim¬ 
ilarly to gg H 7 production in the SM due to charge 
conjugation invariance@ Therefore for this process we con¬ 
sider only scalar, pseudoscalar and axial vector mediators. 
The corresponding total cross-sections are given in Table 
E3 The scalar cross-sections are of the same order but a bit 
larger than the SM HH production [U due to the lack 
of the destructive interference with the triangle with the 
trilinear Higgs coupling. 

The parton-level results for the normalised distribu¬ 
tions of the missing pp for the benchmarks above are 
shown in Fig. l 20 landl 2 Tl for the scalar, pseudoscalar and 
axial-vector mediators respectively. We observe different 
distributions for the scalar scenarios compared to the pseu¬ 
doscalar ones. Indeed the form-factors describing the scat¬ 
tering amplitudes are different depending on the parity of 
the scalar. In the infinite top mass limits, these differ ex¬ 
actly by a factor of 2/3. 



Fig. 20. Missing pr distribution for pp — > XX H for the scalar 
(S) and pseudoscalar (P) mediators. 



Fig. 21. Missing pr distribution for pp — » xxH for the axial- 
vector mediator. 


2 We note here that these selection rules apply only to the 
2 —>• 2 scattering amplitudes. The amplitudes for these pro¬ 
cesses are non-zero when additional QCD radiation is allowed. 


5.3 mono-photon 

For the mono-photon case the only mediator which con¬ 
tributes is the vector one. The contributions of a scalar, 
pseudoscalar and axial-vector vanish due to charge con¬ 
jugation invariance. The cross-sections at 13 TeV for the 
three benchmarks are shown in Table [TJ while the corre¬ 
sponding distributions for the missing transverse energy 
are shown in Fig. [22] Only the resonant scenario gives 
a non-negligible cross-section for unit couplings, implying 
that the mono-photon process in this simplified DM model 
can be accessible experimentally only for limited corners 
of the mass parameter space. 



Pt { “ [GeV] 


Fig. 22. Missing transverse momentum distribution for pp —» 
XXX f° r the three mass benchmarks. The only mediator giving 
non zero results is the vector one. 


6 Conclusions 

We have presented results of the implementation of a sim¬ 
plified dark matter model within MG5 aMC, involving a 
spin -0 or spin -1 mediator, coupling preferentially to the 
top quark. For this model the production of DM particles 
proceeds through gluon fusion loops. We have considered 
the production in association with QCD jets, leading to 
the jets plus missing transverse momentum signature at 
the LHC. The results obtained depend strongly on the 
mass of the mediator and dark matter particles, leading 
to sharper or broader features in the differential distribu¬ 
tions. To provide a reliable description of the distribution 
shapes we have employed matrix-element-parton shower 
merging and matching, and presented results for four me¬ 
diator types: scalar, pseudoscalar, vector and axial-vector. 
Even though we have not considered them here, processes 
with mediators with mixed parity (scalar/pseudoscalar 
and vector/axial-vector) can be also simulated. 

In addition to the jets associated production, we also 
considered the production of DM in association with a Z, 
Higgs and photon. A subset of these production modes are 
forbidden by conservation laws, depending on the media¬ 
tor type. Therefore it is important to complement searches 
with jets with searches for mono-X, as any signal in these 
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Benchmark 

scalar 

pseudoscalar 

axial-vector 

Resonant 

6.98 

i n—2 +34% +1.0% 

‘ 1U -24% -1.2% 

ni o Q +33% +1.0% 
u.ioj _ 23 % _ 12 % 

9 SI in- 2 +36% +1-3% 
z.oi • iu _ 25% _ L4% 

Heavy mediator 

9.31 

in -5 +41% +2.1% 

• 1U _ 2?% _ 2 .i% 

c 7Q in -5 +40% +1.9% 
O.IV • 1U -27% -1.9% 

q ni in-5 +41% +2.1% 
O.U1 • 1U _ 27% _ 2 .i% 

Heavy DM 

1.28 

in -7 +43% +3.0% 

' 1U -28% -2.9% 

2.44 

9 07 1 0“ 5 + 43% + 2 - 9 % 
Z.U/-1U 28% -2.9% 


Table 6. Cross sections (in pb) for gluon induced mono-Higgs production at the LHC at y/s = 13 TeV for the three mass 
benchmarks. A technical cut of 2 GeV has been set on the transverse momentum of all final heavy states but no Higgs branching 
ratios are included. 


Benchmark 

vector 

Resonant 

Heavy mediator 
Heavy DM 

q 1 o in -2+34% +1.1% 

o.io • iu _ 24% _ 12 % 

IQS in— 5 +41% +2.1% 

1.90 • IU _ 27 % -2.1% 

117 1 O — ^ +43% +3.0% 

1.11 • IU 28% -2.9% 


Table 7. Cross sections (in pb) for gluon induced mono-photon production at the LHC at yfs = 13 TeV for the three mass 
benchmarks. A 10 GeV cut is applied on the transverse momentum of the final state photon and + 1 < 2.5. 


channels can straightforwardly provide information on the 
nature of the mediator. 

Our implementation is completely general and public. 
Even though in this study, we presented results at a rather 
inclusive level, the interested reader can use the implemen¬ 
tation to obtain results with the appropriate cuts, mim¬ 
icking the corresponding experimental searches. These re¬ 
sults combined with the experimental measurements can 
be used to extract limits on the various model parame¬ 
ters, also taking into account constraints from direct and 
indirect detection experiments. 
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